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Abstract. The reliable and precise generation of quantum unitary transformations is essential to the 
realization of a number of fundamental objectives, such as quantum control and quantum information 
processing. Prior work has explored the optimal control problem of generating such unitary transfor- 
mations as a surface optimization problem over the quantum control landscape, defined as a metric for 
realizing a desired unitary transformation as a function of the control variables. It was found that under 
the assumption of non-dissipative and controllable dynamics, the landscape topology is trap-free, im- 
plying that any reasonable optimization heuristic should be able to identify globally optimal solutions. 
The present work is a control landscape analysis incorporating specific constraints in the Hamiltonian 
corresponding to certain dynamical symmetries in the underlying physical system. It is found that the 
presence of such symmetries does not destroy the trap-free topology. These findings expand the class 
of quantum dynamical systems on which control problems are intrinsically amenable to solution by 
optimal control. 



1. Introduction 

Central to many problems in quantum control [1] and quantum information processing 2J is the stable 
and precise generation of specific unitary transformations. This task may be viewed as an inverse problem 
where given a desired unitary transformation, one must obtain the values of the control variables of the 
system Hamiltonian whose dynamics generate it. 
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Generally, the connection between the control variables and the unitary evolution operator is suffi- 
ciently complex such that the exact control solution cannot be deduced from first principles. In such 
cases adaptive optimization techniques are commonly applied, such as optimal control theory (OCT) 
methods based on classical variational optimization [3J for computer simulations 4 and optimal control 
experiment (OCE) methods based on evolutionary adaptation [5j for practical laboratory studies p)l. A 
method for the optimal generation of unitary transformations based on OCT methods has been recently 
introduced ^7J, with some promising successful applications in simulation studies [5]. 

A basic question is why such methods have had such a surprising degree of success in such complex, 
high-dimensional problems. A step toward the resolution of this question has been found in recent 
theoretical studies of the quantum control landscape [9l I10| , defined generally as the metric of attainment 
of some optimization objective as a function of the control variables. In the present case, the optimization 
objective is the generation of a specific unitary transformation of the quantum system. 

Consider a controllable quantum system defined on an A^-level Hilbert space whose dynamics are given 
by the time-dependent Hamiltonian Hit). The unitary evolution operator for the system is 



where denotes the time-ordering operation over some time interval [ti,tf]. The optimization prob- 
lem entails choosing the controls within H{t) that steer S{t.i,tf) to some desired target transformation 
W{ti,tf) = W. We assume a degree of system controllability such that for any S{ti,tf) = S selected 
from the unitary Lie group U{N), there exists some choice of controls within H(t) which generates S 
via Eq.(IT]). This full access to U{N) makes it possible to define the control variables of the problem to 
be some parametrization of S itself, in lieu of the controls within H{t). This kinematic representation of 
the control variables is particularly attractive in the sense that the analysis now depends no longer on 
the specific structure of the Hamiltonian, but only on the geometry of U{N). Furthermore, the results 
are generalizable to all finite-dimensional, controllable quantum systems. 

We employ as the distance measure between unitary transformations the squared Hilbert-Schmidt 
metric US' — = 2N — 2KeTT{W'' S), whose nonconstant part we define to be the landscape metric 




(1) 
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function 

J[S]=ReTr{W''S). (2) 

Trace functionals of this type have been studied in a general context [TTl [H] , as well as with specific 
reference to physical [13j and control theoretic applications |14[ I15j . We define the control landscape to 
be the image of the metric function J : U (N) — > R, where each point in U (N) represents some choice of 
controls. We will refer to the space over which the landscape is defined as the landscape domain, which 
in this instance is U{N). 

The optimal control problem of generating unitary transformations of a desired form is essentially a 
search over the control landscape for critical regions at which the first-order variation of the landscape 
function J vanishes. In this sense, the difficulty of the problem is largely determined by the topology of 
the critical points. In a prior analysis of the problem of optimally generating unitary transformations over 
U{N) lOj, it was found that the number of disconnected critical regions of the landscape corresponding 
to distinct values of J scaled as + 1, of which — 1 regions corresponded to local critical points with 
saddlepoint topology, and the remaining to global extrema. Furthermore, it was determined that the 
global extrema solutions comprised zero-dimensional subspaces of U (N) , where the local critical regions 
had the structure of complex Grassmannian submanifolds embedded in U{N) [lOj . 

Although the existence of such local (false) critical points may still act as deleterious attractors for 
optimal searches, the number of such regions grow only linearly with the Hilbert space dimension of 
the underlying physical system. These appear to be competing arguments in the assessment of whether 
this landscape is amenable to adaptive searching methods. In recent numerical studies where genetic 
algorithms, gradient following, and simplex methods have been used, the observed exponential scaling 
of search convergence time suggests that the deleterious infiuences may dominate [IB]. Nevertheless, the 
desirable absence of any local extremum traps over the entire class of such surfaces assures that any 
optimization will eventually succeed in locating a global solution. 

A key open question is whether this desirable trap-free landscape topology is preserved when the 
underlying dynamics of the physical system are restricted or constrained in some manner. In the present 
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analysis, we consider various cases in which specific dynamical symmetries are imposed on the system 
Hamiltonian, specifically relating to total spin, space rotation invariance, and time-reversal invariance. 
In prior analyses [9l llOj , no structure of the Hamiltonian was assumed other than its complex- valued 
Hermiticity (i.e., the imaginary part of the Hamiltonian is nonzero), which corresponds to the symmetry 
class of systems without time- reversal invariance . We presently consider the topology of the control 
landscape defined over physical systems with Hamiltonians from some alternate symmetry classes: 

(1) Systems with time-reversal invariance and integral total spin. 

(2) Systems with time-reversal invariance and space rotation symmetry. 

(3) Systems with with time-reversal invariance, half-integer total spin, without space rotation sym- 
metry. 

The first two classes are described by real symmetric Hamiltonians H ~ Such Hamiltonians arise in 
time-symmetric pulsing strategies in spin control problems [THl UHl 120] > and continuous quantum random 
walks pn [^ . The third class is described by symplectic Hamiltonians H = H^, where the 2N x 2N 
matrix is the symplectic dual 



( .J 




' In ^ 





















(3) 



Such Hamiltonians arise in the dynamics of Gaussian pure states f23J, which have notable applications 
in quantum optics [24 . A more extensive discussion connecting symmetry classes of Hamiltonians with 
their corresponding physical systems can be found in [25[ 126) 

Section 2 presents formal definitions of the landscape function and domain for the case of real sym- 
metric Hamiltonians. From these definitions, the remarkable property of the invariance of the landscape 
topology with respect to the choice of target transformation follows simply. The identification of the crit- 
ical landscape regions as a union of real Grassmannian submanifolds is given in Section 3. The method 
for computing the signature of the critical submanifolds is derived in Section 4. Section 5 recapitulates 
the analysis for the case of symplectic Hamiltonians. Section 6 concludes. 
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2. Landscape Function and Domain 

Consider the class of physical systems whose dynamics are described by real symmetric Hamiltonians. 
Through Eq.([T]), the dynamical propagators generated by such Hamiltonians arc symmetric unitary 
transformations. Let 5 be a programmable N x N symmetric unitary transformation satisfying (i) 
55^ = 5^5* = /jv, where /^r is the N x N identity, and (ii) S = . These constraints leave ^ +^ real 
degrees of freedom in S |26| . For some desired symmetric unitary transformation W, the optimization 
objective is to minimize the landscape metric distance between W and S as defined in Eq.Q. In 
general, optimization heuristics defined over the landscape will seek the critical points of the landscape 
where the first-order variation of the landscape metric function J vanishes. A remarkable quality of this 
optimization problem is that the landscape topology is invariant to the choice of the target transformation 
W. In this sense, the optimization problems for all unitary transformations can be expected to be of 
equivalent difhculty. We presently prove this claim. 

We assume that the range of controls is restricted such that the Hamiltonian generating S is always 
real symmetric. Any symmetric unitary transformation S has the canonical representation S = U'^U 
where U G U{N) 26J. Let Ui,U2 G U{N). The necessary and sufficient condition for = [/J[/2 

is that G 0{N) where 0{N) denotes the real orthogonal Lie group. Therefore, the image of the 

canonical representation mapping U U^U is homeomorphic to the homogeneous space U{N) /0{N) 
[2 7) . Henceforth, we take the space U{N)/0{N) of symmetric unitary transformations to be the land- 
scape domain. 

Since U{N)/0{N) is not closed under multiplication, we must rearrange the argument of the trace 
function J[S] = ReTr (Vw^SVw^), where S and W are symmetric unitary, to ensure that the landscape 
metric is defined strictly over the landscape domain. Noting that Vw^sVw^ ^ 5 is a homeomorphism 
on U{N) /0{N), the image of J[-] = ReTr(-) with the choice of S as the argument is equivalent to 
that with the choice of Vw^SVw^ . Adopting the simpler choice of S as the argument, we obtain the 
topologically equivalent landscape metric function 



=ReTr(5) 
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which has no dependence on the target transformation W . This target-invariant landscape function 
also has the advantage of analytical simplicity, which we adopt for the remainder of the analysis. 



3. Critical Submanifolds 

We presently determine the enumeration and topology of the landscape regions on which is critical. 
We will demonstrate that the number of such regions scales linearly with N and has the structure of real 
Grassmannians embedded in U{N)/0{N). 

Any symmetric unitary transformation can be diagonalized as 



(5) 



where X is an element of the real special orthogonal Lie group SO{N) and is a diagonal operator 



V 



(6) 



of the unimodular eigenvalues {e*'''\ e"^™} of S [28]. Using the cyclic property of the trace, the 
target-invariant landscape metric function simplifies to 



J[S] = ReTr(r2) 



N 



= ^C0S(/3j. 

i=i 



(7) 
(8) 



From Eq.®, we see that the first-order variation SJ^ = — sin(^j dipj vanishes when (pj = ijir, for 

any integers £j. There are + 1 critical values —N,—N + 2, ...,N — 2,N for J', determined by the 
parities of ij,j — 1, TV. The corresponding critical points S = X'^fl^'^'^X comprise equivalence classes 
of transformations orthogonally similar to canonical elements f)'"-' where n denotes the number of even 
integers in the set {£i, ...,£]\[}: 



( 



(9) 



We seek to determine the topology of the iV + 1 equivalence classes of transformations corresponding 
to each distinct critical value. Such equivalence classes are the critical submanifolds of the landscape, 
composed of points at which the first-order optimization condition \7J^ = is satisfied. 

Consider the conjugation of e U{N)/0{N) by some S e SO{N) as a group action G : SO{N) x 
U{N)/0{N) U{N)/0{N) 

g ■ f2(") = s^fi(")s, (10) 

for which SO{N) is the acting group and U{N)/0{N) is the t/-space. The stabilizer (isotropy) subgroup 
STAB(5l(")) of any 57^"^, defined as the subgroup of the acting group SO{N) whose elements map 
back to itself via the ^/-action, is composed exclusively of elements of the form 

/ ^ \ 

'—'71 

S= ^ (11) 

where S„ e SO{n) and S^v-n G SO{N ~ n). Therefore, the stabilizer is a product subgroup of SO{N): 

STAB(r2(")) = SO{n) x SO{N - n). (12) 

Define a mapping which associates a fixed C/-space element $7 with some element of the acting group: 

Gn : SO{N) -> C/(7V)/0(iV), n^Q-^. (13) 

When the domain of is taken to be all of SO{N), the image is simply the orbit of fl: 

0{n) = {E^rtE : E e S0{N)}, (14) 

These orbits, for Q = £7^"^ where n = 0, N, comprise precisely the N + 1 sets of critical points of the 
landscape. 

To obtain their topological structure, we recall that by the orbit-stabilizer theorem, Gq induces a 
bijection 50(A^)/STAB(f7(")) ^ ©(fi^")), which we may sharpen further to be a diffeomorphism because 
SO{N) is a compact Lie group [29^. Therefore, the critical set is composed of the union 

N 

\jGreal{n,N) (15) 
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where 

Greaiin, N) = SO{N)/SO{n) x SO{N - n) (16) 
embedded in U{N)/0{N). The dimensionality of the Grassmannian is well known to be 

N'^-N \n'^-n [N - nf - N + n 



dlU\Greal{n, N) 



(17) 



2 

n{N-n). (18) 



Since the objective of any optimization is to attain the = N ov —N values, corresponding to a perfect 
generation of the desired transformation, an immediate consequence of the foregoing dimensionality 
equation is that the critical submanifolds corresponding to non-global critical have nonzero dimension in 
appropriate subspaces of U{N)/0{N) whereas the critical submanifolds corresponding to global critical 
points, with n = or n = N, strictly have dimension zero. It is therefore of practical importance to 
determine whether the critical submanifolds corresponding to non-global critical points have the topology 
of local maxima or minima, which may act as traps, or of non-trapping saddlepoints. 

4. Hessian Analysis of Critical Points 

A critical point of can be identified as an extremum or a saddlepoint by computing its signature 
(D+, I?-, Dq), an ordered triple of integers denoting the number of upward, downward and flat landscape 
directions at that point. _D+ and Z?_ are commonly referred to as the indices of positive and negative 
inertia, and Dq as the kernel dimension. The Hessian operator of evaluated at S is defined 

n,^ = (19) 

dxi' dxi ' 

where {xi} . n^+n is some set of local coordinates around a point S E U{N)/0{N). Since the 

2—1,..., ^ 

Hessian is symmetric, there exists some coordinate transformations Xi — s- Xi which rotates Ti. into diagonal 
form, where the enumeration of its positive, negative and zero- valued eigenvalues corresponds exactly to 
D_ and Dq. By Sylvester's Law of Inertia, Z3_ and Dq are invariant to changes in the coordinate 
system [50] . 



As a symmetric matrix, Ti is representable as a quadratic form (T\Q{T-L)\T) = ■ QijTiTj, with real 
coefficients {Qij} and {r^jFj} denoting components of a vector |r) of real indeterminates. A corollary 
of Sylvester's Law assures that the coordinate rotation Xi — > Xi induces a transformation of the quadratic 
form into a canonical form of strictly second degree monomials 



(r|Q(H)|r)^^4r? 



(20) 



where the enumeration of the positive, negative and zero-valued real coefficients Qj corresponds to 
D+,D^ and Dq [30] , 

To explicitly connect the Hessian matrix with its quadratic form, let us consider as a mapping 
(i7 o 7) (t) over a parameterized arc 7(t) G U{N) /0{N) satisfying 7(0) ~ S. Taking the second derivative 
of the arc-parameterized Ja.it — Q,we have [3T] 



{Joi)"{Q)={^ 



dxi 



-7^(0) 



dxidxj 



(21) 
(22) 



Noting that ^'^^^^-'^ vanishes since S is critical, we identify the remaining term ^ ^a7^a7^7j('-')7j (0) 
as a quadratic form mapping tangent space vectors 7(0) G Tg U{N)/0{N) to the reals [31]. Transforming 
into the diagonal coordinate system {xi}, we have 



dxj 



(23) 



Associating ^^^^ ^ 



2 7/(0). 

i 

Q -^'i - " Qi and 7^(0) ^ T^, we obtain a direct identification of the Hessian matrix with 



its quadratic form. 

We now compute the Hessian quadratic form (HQF) explicitly. Let us evaluate J" at some critical 
point S in the previously defined representation S = X^ClX. Since S" is a critical point of J^, (l takes 
the form 



/ 



n = 



(24) 



LUN 



where oji = ±l,...,cjAr = ±1. If Cjj = +1(— 1) for all j, we are at a global maximum (minimmn). 
Otherwise, we are at a local critical point. We presently seek to establish that these critical points are 
not local extremum traps. 

To obtain the HQF, we twice differentiate the landscape function argument S along the parameterized 
curve S — V^e*"^* in U{N)/0{N) defined by some real matrix of indeterminates A — A^: 



dt 



Se"^WS = VSiAVS ^ S" ^ -VSA^VS (25) 
t=o 



and evaluate explicitly in terms of the matrix elements of A. 

nj{S) = -ReTi (A^S) (26) 

= -ReTi (xA'^X'^nj (27) 

= -ReTr (^A^qJ (28) 

= -Y,Codl-J2^^Al (29) 

= - Y,(Cj, + - (30) 

where A = XAX"^ and multiplicative constants have been normalized to unity. 

To determine the values of _D_, Dq, we note that for a critical point 5* with eigenvalues +1 and —1 
with multiplicities n and N — n, the only nonzero coefficients in the first summation — X]i<j ('^i + 
of Eg. ([50)1 will be generated by a)i = ujj = +1 or u),; =ujj = —1. There are ways of selecting indices 

i and j satisfying i < j from the n-fold set of indices corresponding to positive values {k : ujk = +1}, 
and thus "^"^^"^ negative-valued monomials Afj in this summation. We add this to the n additional 
negative monomials from the second summation — ^iAf^ to obtain _D_ = " ^" . By similar counting, 
it can be seen that there are £)+ = {N-n)^+N-n pQgj^jyg monomials. Noting that the sum of the indices 
of inertia and the kernel dimension must equal the dimensionality of the domain of the function whose 
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Hessian is being computed 



D++D_ 



+ Do = dimU{N)/0{N) 



(31) 



+ N 
2 



(32) 



it is evident that 



(33) 



D 



2 

+ n 



(34) 



2 



Do 



n{N — n). 



(35) 



Finally, it is evident that for any local critical point for which both N, N — n > 0, both and D_ must 
be nonzero. Therefore, all local critical points are assured to possess saddlepoint topology, and cannot 
act as local extremum traps. 



We consider the class of physical systems whose dynamics are described by symplectic Hamiltonians. 
Since the analysis is almost identical to the real symmetric Hamiltonian case, the discussion will be 
limited to statements of the key assumptions and conclusions. Let S denote programmable symplectic 
unitary dynamical operator generated by such a Hamiltonian via Eq.([T|) with quaternion- real matrix 
elements in the 2x2 Pauli basis {ctq, ax, cry, ^z}- The choice of the Pauli basis for the matrix elements 
induces a 2N x 2N representation of S, which must satisfy (i) SS'' = S^S — I2N and (ii) S = where 
denotes the symplectic dual. Such constraints leave 2N{N — 1) real degrees of freedom for S. 

In analogy to the prior case, we assume that the range of controls is restricted such that the Hamil- 
tonian generating S is always symplectic. Any symplectic dual transformation S has the canonical 
representation S = U"U where U e U{2N) [26j. Let Ui,U2 € U{2N). The necessary and sufficient 
condition for U'[Ui — U2U2 is that t/iC/j"^ G Sp{2N) where Sp{2N) denotes the real symplectic Lie 
group. The image of the mapping giving the canonical representation U — ^ U^U is homeomorphic 



5. Control Landscape Topology for Symplectic Hamiltonians 



11 



to U{2N)/ Sp{2N). Henceforth, we take U{2N) / Sp{2N) to be the space of symplectic unitary trans- 
formations over which we define as the landscape domain for the present case. In cfose anafogy with 
the prior case, we may also define a target-invariant landscape function J^[S] — ReTr(S') defined over 
U{2N)/Sp{2N). 

The first and second order analyses of the critical landscape topology can be replicated for the sym- 
plectic case by using, in analogy to Eq. ^ , the canonical representation for unitary symplectic operators 

S = X^nX, (36) 



with X e Sp(2N) and 



/ 

LUiao 



(37) 



where lui, ...,WAr are the complex, unimodular eigenvalues of S. 

Analogously to the prior case, there are N + 1 critical submanifolds with the structure of symplectic 
Grassmannians. Collectively, the critical set is the union 

N 

[_J Gsympiectic{2n,2N) (38) 

where 

Gsympiecuci2n, 2N) = Spi2N)/Spi2N - 2n) x 5p(2n) (39) 
embedded in U{2N)/Sp{2N), with dimension 

dim G,y,npiecuc{2n, 2N) = N{2N - 1) - (N - n)[2(iV - n) - 1] - n[2n - 1]. (40) 

The HQF can be obtained similarly. Let us evaluate at some critical point S in the representation 
S — X^ilX. In analogy to the prior case, we differentiate S twice over the arc given by 5 ^ 5' = 
= V^iAy/^ S" = -V^A^y/^ defined by some A = (where Aij = a,yCro + Pija^ 



dS 
dt 
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jijCTy +Sijaz with real aij, Pij,"fij,5ij) and evaluate the landscape function with the arc-parameterized 
argument to obtain, up to multiplicative constants 

nj{S) = - ^(^. + co,){a^^ + 4 + 4 + - E (41) 

i<j i 

The number of positive, negative and zero-valued elements of the HQF are, respectively, 

D+ = 2{N -n)^ + 2{N -n) (42) 
^ 2n^ + 2n (43) 
Do = 4iV(iV-n). (44) 
6. Conclusions 

The analysis reveals that the critical submanifolds for the symmetry-restricted landscapes are of two 
types: isolated points corresponding to the global maxima and minima, and Grassmannian submanifolds 
corresponding to the sub-optimal extrema values. Although the sub-optimal Grassmannian solutions 
are more numerous and more voluminous than the global solutions, the Hessian analysis reveals that 
all such local solutions have saddlepoint structure and thus do not act as local traps. Furthermore, the 
invariance of the qualitative landscape structure with respect to the target transformation demonstrated 
in this analysis suggests that no quantum unitary transformation should be any more difhcult to optimally 
construct than another. In addition to prior works [9[|10). the present analysis is a strong endorsement of 
optimal search as a practically viable means of generating unitary transformations for complex systems, 
and demonstrates that the most important qualitative aspects of the landscape are not affected when 
certain symmetry restrictions are placed on the underlying dynamics. 
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